Non-linear transport and heat dissipation in metallic carbon nanotubes 
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We show that the local temperature dependence of thermalized electron and phonon populations 
along metallic carbon nanotubes is the main reason behind this non-linear transport characteristics 
in the high bias regime. Our model that considers optical and zone boundary phonon emission as 
well as absorption by charge carriers is based on the solution of the Boltzmann transport equation 
that assumes a local temperature along the nanotube, determined self-consistently with the heat 
transport equation. By using realistic transport parameters, our results not only reproduce experi- 
mental data for electronic transport, but also provide a coherent interpretation of thermal breakdown 
under electric stress. In particular, electron and phonon thermalization prohibits ballistic transport 
in short nanotubes. 

PACS numbers: 73.63.Fg, 73.23.-b, 65.80.+n 



Carbon nanotubes (CN) are one-dimensional (ID) 
nanostructures that have stimulated broad research inter- 
est because of their unique electrical versatility into semi- 
conductors and metals, depending of their chirality 0. 
From a technological viewpoint, their remarkable elec- 
trical and mechanical properties make them promising 
materials for applications in high performance nanoscale 
electronic and mechanical devices 0, 1^- Among these 
properties, the interrelation between electronic and ther- 
mal transport in these quasi ID structures is particu- 
larly interesting. Early experiments on non-linear trans- 
port in metallic single walled nanotubes (m-SWNTs) us- 
ing low resistance contacts revealed current saturation 
at the 25/iA level, which was attributed to the onset 
of electron backscattering by high energy optical (OP) 
and zone boundary (ZB) phonons in the high bias regime 
. More recently, series of independent high- field trans- 
port measurements on various length m-SWNTs demon- 
strated the absence of current saturation with current 
levels over 60/iA in short samples (^55nm), which was 
interpreted as ballistic transport along the CN |^ |^. 
Among the findings was also the observation of ther- 
mal breakdown and burning under high electric stress. 
In the mean time, the electrical conductance of multi- 
walled nanotubes (MWNTs) under high bias has shown 
step-like decrease caused by the successive burning of the 
CN outer shells In these experiments CNs burn un- 

expectedly at mid-length under stress even on a substrate 
and on the presence of a back gate in a field effect device 
geometry j9i. Despite various attempts to model these 
systems |^ |^ , up to now no coherent interpretation 
has emerged that reconciliates heat dissipation with elec- 
tronic transport and describes thermal effects in m-CNs 
under electric stress. 

In this letter, we show that the non-linear charac- 
teristics of metallic CNs find their origin in the non- 
homogeneous Joule heating along the nanotube, which is 



caused by the thermalized distribution of electrons scat- 
tered by high energy phonons, even in short m-SWNT. 
We specifically show that Joule heating is maximum at 
CN mid-length and, owing to the ID nature of the struc- 
ture, increases drastically with the CN length, resulting 
in thermal breakdown at lower bias than in shorter CNs. 
Our model is based on the Boltzmann transport equa- 
tion with OP and ZB phonon scattering and solved self- 
consistently with the heat transfer equation, providing a 
coherent interpretation of electric and thermal transport 
in m-SWNTs in agreement with experimental data 
In particular we show that the high current level in short 
CNs is not due to ballistic transport but to reduced Joule 
heating. 

We use the linear dispersion relation of electronic states 
around the Fermi level {€{k) — ±hvpk) being vp the 
Fermi velocity. Thus, the Boltzmann equation reads: 
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Here /"(e), F, e and Cph are the distribution function 
for the a energy branch, the electric field, the electron 
charge and the electron-phonon collision integral, respec- 
tively. The index a denotes the energy branches with 
positive (-I-) and negative (-) Fermi velocity in the first 
(1) and second valleys (2) of the m-CN. In metallic sys- 
tems high electron density and strong inter-carrier scat- 
tering thermalizes the electron distribution. We there- 
fore assume that the electron distribution function /"(e) 
obeys Fermi-Dirac statistics with a local electronic tem- 
perature Tei(x): 



/"(e) = 1/(1 + exp((e - e^)/fcsTe,(:r)) 



(2) 



where e^ is the quasi Fermi level of branch a. As a result, 
the collision integral C^^ also depends on the position. 
We neglect acoustic phonons that are only relevant in 
the low bias regime, and consider the contribution of high 
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density OP and ZB phonons 0, ^3 that play a central 
role in energy dissipation in the high bias regime. As 
illustrated in Fig.QIa), the different processes (inter- and 
intra-branch) considered in Cph include both the emission 
and absorption of these phonons with energy (TioJop ~ 
0.2eV) much larger than thermal fluctuations at room 
temperature. Each of these phonons contributes to the 
collision integral as follow: 



C«,(/,r(4 [f\k) (l-rik-q)) -rik) (1 ~nk-q)) 
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where the Greek letter index f3 runs over the two branches 
and two valleys, i stands for OP and ZB phonons and 
{Ra) is the phonon emission (absorption) rate. Hence the 
first (second) two terms in Eq.|21corresponds to processes 
involving the emission (absorption) of a phonon limited 
by Pauli exclusion principle. For instance, the first term 
describes a process in which an electron scatters from a 
state in branch /3 with momentum fc to a state in branch 
a with momentum k — qhy emitting a phonon. In all 
these processes, both total energy and momentum are 
conserved. The emission and absorption rate coefficients 
are given by: 



N 

K{Tl) = 
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Ti exp{huj/kBTL) - 1 



= Ra exp{hu;/kBTL) 



(4) 
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where Tl is the lattice temperature and l/r.; stands for 
the bare scattering rate for OP and ZB phonons that 
we assume to be independent of carrier energy in a first 
approximation. In computing the collision integral, we 
make the key "ansatz" that electrons and lattice are in 
local thermal equilibrium (i.e. TL — TL{x)—Tei{x)). We 
define the electron density as: 
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r{k)dk 



(6) 



where Ec is the bottom of the conduction band, and the 
current to be: 



(7) 



where the +(— ) index corresponds to the branches with 
positive (negative) Fermi velocity. Then integrating Eq.^ 
over the momentum for each branch, and properly ac- 
counting for different branches, we obtain: 



vpdx (ri^ — n ) = 



(8) 
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VFd^{n++rr)~^ =2 jdkCphil ,T{x))^ 2Cph{I,T{x)) (9) 



Eq. ISjis the expression of the current conservation in the 
system which by symmetry e^.^ = ^-nd with charge 
neutrality in the CN yields e^^^ 2) — 2)- Integrat- 

ing Eq.|n|over the length of the nanotube L, and assum- 
ing equal electron densities at the contacts, we find the 
voltage drop Vds along the nanotube given by: 

Vds = / dxCph{I.T{x)) (10) 

This equation implicitly depends on the current and the 
temperature profile along the nanotube, and must be 
solved self-consistently with the heat transport equation 
to obtain both current and temperature profile. 
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FIG. 1 
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(a) Scattering processes considered in the calcu- 



The intravalley (left) and intervalley (right) tran- 
sitions with emission and absorption of optical phonons are 
included, (b) IV characteristics for different constant tem- 
peratures along the tube. Dashed line: asymptotic behavior 
for all temperatures in the high bias regime. 

As a particular case, it is interesting to compute the 
IV relation from Eq. IIUI by assuming only OP phonon 
scattering at constant temperature (i.e. T{x) — Tq) in the 
CN. The results are plotted in Fig. where the high 
bias regime exhibits an asymptotic behavior independent 
of temperature which is given by: 



Vds{I) 



1 L 

Go vpT, 



(11) 



op 



where Ii^^ =euJop/Tr is the threshold current, correspond- 
ing to the onset of electron backscattering by OP phonons 
Q and Go = 2e^//i is the quantum conductance. We 
point out that Eq. 1111 is not consistent with the cur- 
rent interpretation of electrons accelerated ballistically 
in the electric field until acquiring enough energy to emit 
a phonon, but rather results from the imbalance between 
the population of the energy branches with positive and 
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negative Fermi velocities. By considering ZB phonons, 
the voltage drop in the m-SWNT is a linear combination 
of expressions similar to Eq. IllL which would still result 
in a threshold current, but with a more complicated ex- 
pression. 




/ [^lA] 

FIG. 2; Temperature difference between the middle of the 
tube and the leads as a function of the current for different 
CN lengths 5|. Inset: Temperature profile along the NT for 
different values of FL. 

When heating effects become relevant, thermal dissi- 
pation is taken into account self-consistently with Eg. 1101 
We consider that two mechanisms for the heat dissipa- 
tion: i) diffusion through the supporting substrate (if the 
CN stands on one), and ii) flow through the contacts. 
Hence, defining AT = T{x)—Tq (where Tq is the tem- 
perature of the substrate and leads), the heat equation 
becomes 1131: 



(12) 



where k is the thermal conductivity, 7 is the coupling co- 
efficient with the substrate and q* is the power dissipated 
per unit volume. Here we make the usual approximation 
that process (i) is proportional to the local temperature 
difference between CN and substrate. In our calculations 
both the thermal conductivity and the coupling coeffi- 
cient are assumed to remain constant along the tube. 
The coefficient 7 is given by: 
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where Ksu6, t and d are the thermal conductivity of the 
substrate, the diameter of the nanotube and the thickness 
of substrate, respectively. We also assume that the power 
is homogeneously generated along the CN and given by 
Joule's law: 

q*=jF (14) 

where j = I/A is current density through the effective 
cross section A and the electric field F is given hy F — 
\yDs/L\- Then the solution for temperature profile is 
given by: 

cosh(Fa;) 



AT{x) 



1 - 



cosh(FL/2) 



where F = j n. 

Different scenarios can take place depending on the 
value FL (see the inset of Fig. EJ. On the one hand, 
for FL <C 1 diffusion through the substrate is negligible 
and the temperature profile exhibits a parabolic shape. 
On the other hand, for FL^l, heat basically dissipates 
through the substrate and the temperature is almost con- 
stant along the CN. This latter situation occurs in long 
tubes strongly coupled to the substrate. In any case, the 
highest temperature point is at the middle of the tube. 
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FIG. 3; Comparison between theoretical and experimen- 
tal IV characteristics for different CN lengths: (a) Ref. 
(b) Ref. Q. 

In our calculations we use To=300K. The energies of 
OP and ZB phonons are ?ia;op=0.20eV and ?ia;^;,=0.16eV, 
respectively. We use the standard accepted value for 
the thermal conductivity k (30W/cmK) and F = 
10^^Wcm~'^K~^. Fig. 121 shows the temperature differ- 
ence in the middle of the tube (AT(0)) as a function of 
/ for different CN length, corresponding to the data of 
Ref. j^. The longer the nanotube, the faster the rise in 
temperature as the threshold current is overcome. This 
is due to the fact that dissipation occurs over a longer 
distance while heat removal mainly takes place at the 
contacts in ID structures. Estimates for the breakdown 
temperature correspond to 800° C Therefore, short 
tubes are expected to carry larger currents before ther- 
mal breakdown. As shown in Fig. O the results for the 
IV characteristics are in good agreement with the exper- 
imental data 1^ 1^ . Deviations in the low bias regime are 
mainly due to the absence of acoustic phonons scattering 
in our model. For the sake of simplicity, we assume re- 
laxation times for OP and ZB phonons with equal values, 
which are t= (13±2)fs for the first and r = (6.9±1.5)fs 
for the second 0| sets of experimental data. The differ- 
ence between the two values could be due to the fact that 
CNs may have different diameters with different phonons 
spectra (breathing modes) in each case, and contact 
quality. Nevertheless, the obtained mean free paths (be- 
tween 6nm and lOnm) are consistent with experimental 
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previous estimates0, @ . For long tubes (i.e. L=700nm 
in the first case |5| and for L >1000nm in the second 
case 0), we need to increase in both cases the relaxation 
times to (27±3)fs to fit the experimental data. Since 
dissipation is considerably stronger in long CNs, these 
longer times may be associated to the emergence of non- 
linear thermal effects in m-SWNT thermal conductivity 
not taken into account in Eq. 1121 and for which the tem- 
perature and geometry dependence at room temperature 
or higher remains an open issue [rRll8j|. 
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FIG. 4: IV characteristics for different tube lengths and 
diameters. Short tubes (lOOnm and 300nm) are practically 
independent of the diameter, while in the long tube (700nm), 
deviations appear for the small diameter (Inm). 

Finally, Fig. 0] shows the /1^-characteristics obtained 
for CNs of different diameters and lengths, assuming that 
both the thermal conductivity and the relaxation time 
are equal among tubes. Despite this strong assumption, 
the relevant issue to emphasize here is the weak depen- 
dence of the IV characteristics on the size of the CN. Ap- 
preciable deviations can be observed only in 700nm tube 
(dashed line) for small diameters (~lnm). This weak 
relation is consistent with the interpretation of Collins 
et al. 0, 01 for the electrical breakdown under electrical 
stress observed in MWNTs, where different layers in a 
MWNT (separated by about 0.4nm) carry similar cur- 
rents in the high bias regime. The breakdown of suc- 
cessive carbon layers produces approximately constant 
diminutions of the current in the high bias regime because 
/V^ characteristics are geometry independent. Moreover, 
the highest temperature arises at the CN mid-length (in- 
set of Fig. 13) and therefore, electrical breakdown is, as ex- 
perimentally observed, expected to take place there too. 

In conclusions, we have shown that the consideration 
of a thermalized electron distribution in local equilibrium 
with a non-homogenously heated lattice through OP and 
ZB scattering determined self-consistently by the cur- 
rent level account for the non-linear IV characteristics 



of the m-SWNTs in the high bias regime. The magni- 
tude of the temperature variation as a function of the 
CN lengths is consistent with the occurrence of thermal 
breakdown at mid-length for long CN under electrical 
stress. While the dependence of thermal conductivity on 
temperature still reamins under investigation, our self- 
consistent model provides a coherent picture of the onset 
of thermal effects with electronic transport in m-SWNT. 

The authors are indebted to A. Cangellaris for fruitful 
discussion. This work was supported by the Beckman 
Institute for Advance Science and Technology and NSF 
- Network of Computational Nanotechnology. 
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